Relativistic Attosecond Electron Bunches from Laser-Illuminated Droplets 
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The generation of relativistic attosecond electron bunches is observed in three-dimensional, rela- 
tivistic particle-in-cell simulations of the interaction of intense laser light with droplets. The electron 
bunches are emitted under certain angles which depend on the ratios of droplet radius to wavelength 
and plasma frequency to laser frequency. The mechanism behind the multi-MeV attosecond electron 
bunch generation is investigated using Mie theory. It is shown that the angular distribution and the 
high electron energies are due to a parameter-sensitive, time-dependent local field enhancement at 
the droplet surface. 

PACS numbers: 52.38.Kd, 42.25.Fx, 52.65. Rr, 52.27.Ny 



Oh. 



in 
O 

Oh. 



> 



O 

o 



The understanding of the laser energy conversion into 
fast electrons and ions is of utmost importance for the de- 
sign of efficient "table-top" particle accelerators [H, 0, [1] . 
The laser plasma-based acceleration schemes discussed 
so far in the pertinent literature may be divided into 
two groups according to whether the plasma is under- 
dense, i.e., the plasma frequency tOp is smaller than the 
laser frequency uj, or vice versa. Wake- field accelerators 
and the so-called "bubble-regime" (see, e.g., 0, H, 0|) 
fall into the former category while the interaction of in- 
tense laser pulses with solid surfaces or thin foils (see e.g., 
0, i, B E El) belongs to the latter, overdense regime. 

Of particular interest are finite-size targets where fast 
particles cannot escape into the field-free bulk material 
but yet the density of the accelerated particles may be 
sizable. In recent years the nonrelativistic interaction of 
intense laser light with small, subwavelength-size clus- 
ters has been thoroughly investigated In particular, 
the efRcient absorption of laser energy, leading to high 
charge states and thus to intense line emission as well 
as to energetic electrons and ions was studied. In such 
small clusters of radii R < 6 <^ X, where A is the laser 
wavelength and i5 is the skin depth of the cluster plasma, 
the effect of the cluster on the propagation of the laser 
pulse needs not be taken into account. Although the 
plasma, which is created via field ionization on a sub- 
cycle time scale, is overdense, screening of the cluster 
interior only occurs due to polarization but not due to 
a skin effect. Technically speaking, the dipole approxi- 
mation can be applied to the nonrelativistically intense 
laser field, E{r, t) ~ E{t). As a consequence the electron 
dynamics mainly occurs in the laser polarization direc- 
tion while the v x B-force in laser propagation direction 
and the influence of the scattered electromagnetic field 
on the particle dynamics can be safely neglected. 

The other extreme of intense laser-matter interaction 
is constituted by targets of sizes much larger than a wave- 
length, e.g., a laser beam impinging on a solid surface. 
The absorption mechanisms in this case have also been 
extensively investigated [l3| . It is well-known, for in- 




FIG. 1: (Color online) Electron isocontour surfaces (1% of 
rieo) of a _R = A/4 Ife-droplet of density n^o = 22 Uc with = 



in a 16-cycle sin -laser pulse (A — 800 nm) of 
W/cm^ (a = 1) at t = 10 cycles. The laser 



cm 

intensity 2 x 10 
pulse propagates in z-direction, the electric field is directed 
along X. 



stance, that for perpendicular incidence electrons are ac- 
celerated in laser propagation direction (i.e., into the bulk 
material). 

In the present Letter we will focus on electron accel- 
eration in the regime where R and A are of the same 
order of magnitude (i.e., rather droplets than clusters) 
and the laser intensity is relativistic, i.e., the ponderomo- 
tive energy Up exceeds the electron rest energy mc^ . By 
changing the droplet radius from i? <C A (small-cluster 
limit) to i? A (solid-surface limit) the emission angle 
of electrons is expected to decrease from 6* = 7r/2 to 
with respect to the propagation direction k. Moreover, 
the emission is expected to occur in the [E, fe)-plane for 
a linearly polarized laser field of amplitude E. 

For a sphere in a plane electromagnetic wave all an- 
gles of incidence occur simultaneously. It is known from 
laser-plasma interaction studies that, depending on the 
plasma scale length and the laser intensity, there exists 
an optimal angle of incidence for the absorption of laser 
energy the steeper the plasma gradient the more 

this optimal angle is pushed towards 7r/2, i.e. grazing 
incidence. For sufficiently large scale lengths resonance 
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absorption (see, e.g., [l[) can occur most efficiently under 
a certain optimal angle. Indeed, this effect, resulting in 
the electron emission opposite to the incoming laser pulse 
has been observed in experiments with droplets and two- 
dimensional model simulations [15|. 

For R > 6 the self-consistent electromagnetic field 
needs to be calculated. Numerically we do this by 
means of three-dimensional, electromagnetic, relativis- 
tic particle-in-cell (PIC) simulations. In such simu- 
lations the particle dynamics, the deformation of the 
target, and the corresponding modification of the field 
is automatically taken into account. Analytically, the 
self-consistent electromagnetic field around and inside a 
laser-illuminated droplet of given dielectric permittivity 
e can be calculated using the corresponding solutions 
to Maxwell's equations put forward by Clebsch, Lorenz, 
Mie, Debye, and others. The electromagnetic scattering 
by dielectric or metallic spheres is commonly called "Mie 
theory" ITtI . Of course, the self-consistent par- 
ticle dynamics far from equilibrium in very intense laser 
fields cannot be captured by a dielectric constant and 
thus is not included in Mie theory while it is included 
in the PIC simulations. However, we will show that Mie 
theory is nevertheless capable of explaining the angular 
distributions observed in the PIC simulations and the 
high electron energies exceeding the expected few times 
the ponderomotive energy. 

We start by presenting typical results from PIC simu- 
lations of the interaction of an intense plane wave laser 
pulse with a pre-ionized He droplet. Figure [T] shows a 
snapshot of the electron isocontour surfaces correspond- 
ing to 1% of the initial electron droplet density. The 
dimensionless vector potential amphtude a = \eA/mc\ — 
\eE/mujc\ of the 16-cycle sin^-laser pulse (A = 800 nm) 
was a = 1, corresponding to a laser intensity / ~ 
2 X 10^^ W/cm^. Electron bunches emitted each half cy- 
cle under plus/minus a certain angle 9 (with respect to 
the propagation axis z) in forward direction are clearly 
visible. The bunches are mainly confined to the plane of 
incidence (i.e., the (a;, z)-plane). 

The electron energy and density for the same He- 
droplet in a four-times more intense laser pulse is pre- 
sented in Fig. [21 Figure [2^ shows that the bunches 
consist of electrons with energies up to ~ 6MeV. The 
ponderomotive energy Up = mc^{^Jl -\- a? /2 — 1) is only 
~ 0.37 MeV. As the spatial width of the bunches is much 
smaller than a wavelength, the temporal bunch structure 
is already in the attosecond domain. Similar bunches 
have been observed in PIC-simulations of the interac- 
tion of ultraintense laser pulses with solid surfaces under 
grazing incidence or with plasma channels [l9| . Such 
bunches have potential applications in attosecond elec- 
tron diffraction experiments or the generation of coher- 
ent short- wavelength radiation via Thomson scattering. 

In the following we will explain the multi-MeV electron 
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FIG. 2: (Color online) Kinetic electron energy (a) and density 
(b) contour plots (in the y = 0-plane) of the He-droplet of 
Fig. [T] in a 16-cycle sin'^-laser pulse (A = 800 nm) of intensity 
8 X 10^*^ W/cm^ (a = 2) at f = 8 cycles. 

energies and their emission angles using Mie theory. The 
latter gives the analytical electromagnetic field config- 
uration around and inside the (unperturbed) droplet in 
terms of series expansions over Legendre and Bessel func- 
tions [i6|] . In the case of interest to us the droplet is con- 
ducting (e = 1 — cjp/cj^), the skin depth S ~ c/iVp <^ R, 
and the electric field on the droplet surface is perpendicu- 
lar to it and quickly decays inside the droplet. Electrons 
at the droplet surface can be pulled out of the droplet 
or pushed inside, depending on the time they appear at 
the surface [3] . Figure [3] shows the absolute value of the 
radial electric field for an ui^juP' = 22 times overdense 
droplet at the droplet surface vs angle 9 and time for 
R — A/4 and R = A/2 according Mie theory. Note that 
in our Mie-calculations the incident plane wave field is of 
the form Einc,xiz,t) = E cos{kz — ojt), i.e., has a constant 
amplitude. Figure [3] shows that Mie theory predicts cer- 
tain times and angles at which the electric field at the 
surface is largest. Moreover, these electric field maxima 
may exceed the field amplitude E (by a factor 3 in Fig.[3K 
and 2.3 in Fig. [3]d). For the bigger droplet the optimal 
angle is smaller, which means that the maximum electric 
field occurs more in propagation direction 9 = Q. Sec- 
ondary local maxima are visible at ~ 0.457r. 

The two insets in Fig. [3] show the electric vector field 
in the y — 0-plane according to Mie theory at the "opti- 
mal" times when the electric field at the surface is highest 
and points inwards (i.e., it pulls electrons outwards). In 
Fig. [3^, for instance, t ~ 0.275 eye. and the largest field 
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FIG. 3: (Color) Absolute value of the radial electric field for 
a 22 times overdense droplet at the droplet surface vs angle 
6 and time, as predicted by Mie theory for R — X/4: (a) and 
R = A/2 (b) in the y — 0-plane. The color indicates the 
electric field in units of the incident field strength E. The 
insets show the electric vector field in the y = 0-plane at 
the "optimal" times (white arrows) when the electric field 
at the droplet surface is highest and pointing inwards, i.e., 
t — 0.275 eye. (a) and t = 0.675 eye. (b), respectively. 



occurs under the angle 6 = 0.277r. It is clear that under 
the same angle but x < the electric field has the same 
absolute value but points outward. As a consequence, 
electrons entering from inside the droplet into such a field 
configuration under this optimal angle will be efficiently 
accelerated outwards (inwards for x < 0). Half a laser 
cycle later the situation reverses and electrons at x < 
will be accelerated outwards (inwards for a; > 0). This 
explains the observation of alternating electron accelera- 
tion into two directions each laser cycle. 

Figure H] shows angle-resolved ion and electron energy 
spectra obtained from the PIC simulations. The spectra 
were taken at time t = 9 cycles when the fast electrons 
are still inside the simulation box. The ions are hardly set 
into motion at such an early time. However, the angular 
distribution of the ions in Figs.S^ andSj; have imprinted 
on them already the anisotropy due to the field distribu- 
tion at the droplet surface. In fact, the emission angles 
are more easily inferred from the ion distributions than 
from the electron distributions in Figs.SJj andlHi because 
the electrons change their direction as they move away 
from the droplet. This is the reason why in Figs.[3)3 and 




-0,03 0,00 0,03 -■^ -i 1 2 

E., MeV Ee, MeV 



FIG. 4: (Color online) Angle-resolved ion and electron kinetic 
energy spectra for the R — X/i droplet [panel (a) for ions and 
(b) for electrons] and the R — X/2 droplet [panel (c) for ions 
and (d) for electrons] after t = 9 cycles in a 16-cycle sin^-laser 
pulse (A = 800 nm) of intensity 2 x lO" W/cm^ (« = !)■ The 
emission angle 9 [indicated in panel (c)] was determined as6 = 
arctan(pa;/pz) with p the momentum. The kinetic energy for 

electrons and ions reads iSe.i = [m^ iC'^+p'^c^]^^^ —m^AC^ ■ The 
color-coding is proportional to the logarithm of the particle 
number. 

|4ti the energetic electrons are aligned stronger in forward 
direction than the slow electrons. 

The emission angles inferred from the PIC results in 
Fig. H are e = 0.277r for the R ^ A/4-droplet and 
6 = 0.147r for the R = A/2-droplet with a secondary 
maximum at OAn. These results are in excellent agree- 
ment with the Mie results presented in Fig. \3\ We have 
performed a systematic study, comparing the emission 
angles predicted by Mie theory with those from the PIC 
calculations (a = 2) for droplet radii between A/8 and 
A/2, and in all cases found very good agreement. 

The unexpectedly high electron energies can be ex- 
plained by the field enhancements at the droplet sur- 
face. As shown in Fig. [3^ the field amplitude at the 
surface of the R = A/4-droplet is 3 times the incident 
field amplitude. We have performed test-particle calcu- 
lations where we placed electrons with zero initial veloc- 
ity on the droplet surface at various phases with respect 
to the analytical, time-dependent Mie field configuration. 
The maximum kinetic energies acquired by such test elec- 
trons are in good agreement with the maximum energies 
observed in the PIC simulations. 

In the limit of small plasma droplets kR = 2ttR/X <^ 1 
the maximum radial field at the droplet surface predicted 
by Mie theory is 

ujp/ijj^ — 3 
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FIG. 5: (Color online) Same as in Fig. [2] but for a = 20 
(8 X 10^*^ W/cm^). The arrow in panel (b) indicates an elec- 
tron bunch which, after the pulse, acquired an energy of 
~ 130 MeV. 



For iOp > one obtains eI"^'^^^ / E = 3, i.e. a three- 
fold field enhancement, as observed in Fig. [3^ (although 
kR = 7r/2 > 1 in this case). For, e.g., ujp = 2uj, the pre- 
dicted field enhancement according H]) is eI"^'^^^ /E — 9 
so that one could argue that even more energetic elec- 
tron bunches can be produced for lower-density droplets. 
However, only for a <C 1 do the PIC calculations repro- 
duce the predictions of Mie theory as far as the field en- 
hancements at the droplet surface are concerned. This is 
because for relativistic laser intensities an only few-times 
overdense droplet quickly dissolves during the rising edge 
of the laser pulse. 

Before concluding we show the simulation result for the 
droplet of Fig. [2] in an ultra-intense laser field of intensity 
8 X 10^°W/cm^. At the plotting time {t = Scycles, i.e., 
at the maximum of the pulse) the highest electron energy 
observed is ~ 50 MeV. However, if one follows the elec- 
tron bunch indicated by an arrow in Fig. [5Jd in time one 
finds a final energy of ~ 130 MeV, i.e., 20 times the pon- 
deromotive energy. From Fig. one infers that almost 
all electrons are removed from the droplet. Nevertheless 
the angle under which the electron bunches are emitted 
is well-described by Mie-theory. However, the energetic 
electrons are not only located inside the electron bunches 
anymore at such high intensities, as is clearly visible in 
Fig. [5ji. 

In summary, we showed that multi-MeV attosecond 
electron bunches are produced when intense laser fields 
interact with overdense droplets of diameters comparable 



to the laser wavelength. The attosecond electron bunches 
are emitted each half laser cycle under plus/minus a cer- 
tain angle in the polarization plane. The preferred elec- 
tron emission angles and the high kinetic energies arise 
due to local field enhancements at the droplet surface 
that can be calculated using Mie theory. Relativistic 
attosecond electron bunches may be used for the gen- 
eration of short-wavelength radiation via scattering of a 
counter-propagating laser pulse, time-resolved structural 
imaging, or plasma diagnostics. 
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